Non-uniform chiral phase in effective chiral quark models 
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^ ■ We analyze the phase diagram in effective chiral quark models (the Nambu- 

Jona-Lasinio model, the cr-model with quarks) and show that at the mean-field 
level a phase with a periodically-modulated chiral fields separates the usual phases 

■ with broken and restored chiral symmetry. A possible signal of such a phase is the 
production of multipion jets travelling in opposite directions, with individual pions 
having momenta of the order of several hundred MeV. This signal can be interpreted 

■ in terms of disoriented chiral condensates. 
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The problem of phase transitions in hadronic matter is one of the most impor- 
Q-h! tant issues in strong-interaction physics. It has been studied theoretically for 

many years with a variety of methods, with the results gaining importance at 
the approach of the operation of RHIC. We now believe that in the (T-fi) dia- 
gram, where T denotes the temperature and \i the baryon chemical potential, 
one can recognize a quite complicated structure: phases with broken/restored 
chiral symmetry, confined/deconfined phases, color superconducting phase (in 
the case of two light flavors) PJ, or color-flavor locked phase (in the case of 
three light flavors) p|. In certain limits detailed features of phase transitions 
are known from first-principle QCD calculations (color superconducting gaps 
at asymptotically- large baryon densities 0), or general symmetry considera- 
tions (e.g. the universality class of the chiral phase transition M in the massless 
limit). Lattice calculations provide information for finite-temperature QCD at 
zero baryon density fl5|. However, these "exact" methods are inapplicable to 
the case of moderate baryon densities (up to a few times the nuclear satura- 
tion density), where out of necessity we have to rely on models In fact, 
a lot of our experience comes from model calculations, just to mention the 
expectation that chiral symmetry is eventually restored as the baryon density 
is increased. 



In this letter we consider the chiral phase structure for a class of low-energy 
chiral models: the Nambu -Jona-Lasinio model || and the linear cr-model with 



quarks 0. These models have been used extensively over the past years to 
describe low-energy properties of mesons and baryons. They have also been 
used to describe dense and hot system of the Fermi gas of quarks [[]], in 
particular to study chiral restoration at high T or fi. What is usually claimed 
is that at the mean-field level (synonymous to one-quark-loop or leading- iV c ) 
the model has two phases: with broken chiral symmetry, i.e. large dynamically- 
generated quark mass M, which appears at low T and [/,, and with restored 
chiral symmetry, with massless (or almost massless) quarks, at higher values 
of T or jj,. We argue that the situation is much more complicated if one allows 
for more general mean-field solutions of the model. We show that for a certain 
range of T and [i there appears a non-uniform phase with broken rotational 
and isospin symmetry. In the T-\x diagram this phase separates the usual 
chiral-symmetry-broken phase from the chiral-symmetry-restored phase. The 
results for the linear a-model with quarks at T = were obtained a long time 
ago by Kutschera, Kotlorz, and one of us (WB) [1(J. Here we generalize the 



calculation to finite temperatures, as well as apply it to the Nambu-Jona- 
Lasinio model. 

The Lagrangian densities of the a-model and of the partially-bosonized Nambu- 
Jona-Lasinio model have the form r/> [ifft — g (a + i'j^T • n)] tp + L m (a, tt), where 
ip denotes the quark fields, g is the quark-meson coupling constant, and the 
mesonic part L m (a, 71") is specific to the model. We consider spatially non- 
uniform, time-independent chiral fields: 

M M 
o(x) = — cos(g ■ x), tt (x) = — sin(g ■ x), ir±(x) = 0, (1) 

9 9 

where M and q are parameters, which are determined dynamically. Such an 
ansatz has been considered for the first time in Refs. JTTJ in context of the pion 
condensation in nuclear matter. In the case of vanishing q we recover a = 
7? = 0. The Dirac equation with fields (0) can be solved exactly HTT| , |T0| , which 
facilitates greatly practical fliO| , |P2| and theoretical |T3] applications. The quark 
spectrum has two branches: 



E ± (p) = ^M 2 +p 2 + -q2±^/M 2 qi + (q-p) 2 } (2) 

where p labels the momentum of the quark. The basic quantity of our study 
is the grand thermodynamical potential density fi(T, A 4 ! M, q). Its global min- 
imum with respect to M and q at fixed T and fi determines the ground state 
of the system subject to conservation of baryon number. We stress that our 
calculation is self-consistent, i.e. both quark and meson fields are solutions 
to the equations of motion. At the mean-field (one-quark-loop) level we have 
explicitly 
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n(T,n;M,q)=n (M,q) (3) 
d 3 p 



2iV c Tg/^l„[(l + 



e (ti-Ei(p))/T\ f x + e -^ +Ei {p))/T 



The medium part of Q describes the ideal gas of quarks with the spectrum 
(0). The vacuum part, f2 , depends on the model considered. In the cr-model 

m 2 F 2 , x 2 F 2 , s 

fiJ(M, q) = (M 2 /M 2 - l) + -f (M 2 /M 2 ) g 2 , (4) 

where the parameter M = gF n is the vacuum value of the constituent quark 
mass, F n = 93MeV is the pion decay constant, and m a is the vacuum value 
of cr-meson mass, treated as a model parameter. The first term in (Q) comes 
from the Mexican Hat potential, while the second term is the kinetic energy 
of the meson fields. In the Nambu-Jona-Lasinio model we have 

Q™ h (M, q) = (M 2 - M 2 ) /{AG) + V(M) + 5 (m 2 /M 2 ) f + 0(q 4 ), (5) 

where G is the four-quark coupling constant, and V(M) is the Dirac sea contri- 
bution to the energy for the case q = 0. Here we adopt the simple 3-momentum 



regulator, hence V(M) = -2N c /it 2 J a dpp 2 \^/p 2 + M 2 - ^p 2 + M 2 J , where 

A is the cut-off parameter. In the calculation presented below we use the pa- 
rameters of Ref. @: G = 5.01GeV" 2 , A = 650MeV. Equation (|) has been 



gradient-expanded in powers of q. The coefficient of the quadratic term is 
model independent. Higher-order terms depend on the value of A. However, 
their contribution is small for moderate values of q [H| and we can safely 



drop them. With this simplification the two models become very similar. The 
same conclusion holds for other variants of the Nambu-Jona-Lasinio model, 
which differ by the regularization method. In fact, the basic difference between 
the models is the difference of values of Qq at the minimum and at the point 
M = 0, q = 0. The larger this difference, the more difficult it is to restore 
chiral symmetry. 

The essence of the dynamics of our system can be understood as follows: as we 
increase the chemical potential, it becomes favorable for the system to develop 
a non-zero q. This is because the quark energies on the E_ branch lower with q, 
and the energy gain overcomes the repulsion of the meson kinetic term in Qq. 
Thus, for certain T and /i the ground state has finite q. At the further increase 
of fi the mass of the quark drops to 0, and chiral symmetry is restored. The 
phase diagram for the Nambu-Jona-Lasinio model with parameters of Ref. 
TA\ is shown in Fig. 1(a). We can see three phases: at low fi and T we have 



the uniform chirally-broken phase, with M > and q = 0, wrapped around it 
is the phase with M > and q ^ , and at high \x or T we find the chirally- 
restored phase, with M = 0. The lower line in Fig. 1(a) describes the first-order 
phase transition between the uniform and non-uniform chirally broken phases. 
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Fig. 1. (a) The T-fj, phase diagram for the Nambu-Jona-Lasinio model with pa- 
rameters of Ref. [|i~4}| . The lower line separates the uniform and non- uniform chi- 
rally-broken phases (first-order phase transition) , and the upper curve separates the 
non-uniform chirally-broken phase from the chirally-restored phase (second-order 
phase transition), (b) Same as (a) in the T-baryon-density space. The ares between 
the two lower curves is the region of equillibrium of the uniform and non-uniform 
chirally-broken phases, (c) Values of \q\ and M in the non-uniform chirally-broken 
phase at the phase transition (dotted and dashed lines, respectively), and the value 
of M in the uniform chirally-broken phase at the phase transition (solid line). 

The discontinuities in the M and q parameters can be read-off from Fig. 1(c). 
Correspondingly, in the baryon density - temperature diagram (Fig. 1(b)) 
we notice the region of equilibrium of the two phases (the region between 
the two lower lines). The upper curve in Fig. 1(a) shows the second- order 
phase transition between the non-uniform chirally broken phase to the chiraliy- 
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restored phase. At this phase transition the order parameter M vanishes at 
non-zero q. Interestingly, the character of the transition between the uniform 
and non-uniform chiral phases depends on the parameters of the models. In 
the cr-model for m a = 800MeV we find a very similar diagram as in Fig. 1. 
However for m a = 1200MeV the phase transition is second-order. It changes 
character around m a = 1040MeV. 



There is always the question of applicability and reliability of a model cal- 
culation like ours. Firstly, we need to have sufficiently high baryon densities 
or temperatures to trust the calculation. Physical systems at low densities 
consist of nucleons, and we need to dissolve them into quarks for the model 
to be applicable. By geometrical arguments we may expect this to happen at 
densities of the order of a few nuclear-matter saturation densities. Note that in 
Fig. 1 this is the region of the non-uniform chirally-broken phase, thus we may 
hope that the model is valid there. Secondly, the size of the system should be 
sufficiently large such that our infinite-size calculation is meaningful. It should 
be larger than 2n/\q\ in order to accommodate at least one period of the chiral 
field. Using the values of \q\ ~ 200 — 600 MeV this gives the minimum size 
in the range 2 — 6 fm, easily accessible in heavy-ion collisions. In addition, 
a more elaborate calculation should incorporate meson fluctuations. At low 
temperatures pions are easily excited due to their small mass, and such effects 
should certainly be included. Note, however, that inclusion of meson loops in 
the study of the quark condensate in Ref. ]nj modified, but did not invalidate 
the mean-field results. Finally, the effect of the explicit breaking of chiral sym- 
metry by the finite current quark mass, m, should be incorporated. In studies 
with q = H this effect was not quantitatively small. Moreover, with finite 
m chiral restoration is no longer a second-order phase transition, but becomes 
a smooth cross-over. Another question is whether within the mean-field ap- 
proximation the solution with fields (|l|) is the ground state. What we have 
verified is that for a certain range of T and fi we have a lower-energy state 
than the uniform phase. A priori, there may exist yet lower-energy, up-to-now 
unknown, mean-field solutions. 



We conclude with a discussion of a possible signal of the non-uniform chi- 
rally broken phase. Suppose that sufficiently large domains of the non-uniform 
chirally-broken phase are created in the cooling process of the plasma formed 
in a heavy-ion collision. Such a domain can be treated as a classical source 
of pions, in full analogy to disoriented chiral condensates (see e.g. |16|). The 



coherent decay of such a source would lead to a characteristic signal in the 
form of two multi-pion "jets" , with pion momenta peaked around ±g (in the 
rest frame of the source), where q is taken at the phase-transition point, and 
is typically 200 — 300MeV. Indeed, according to Bjorken's description, pion 
quanta 5(p a (x) satisfy the equation (□ + m^) 5(p a (x) = j a (x), where j a (x) is 
the source, in our case proportional to the classical pion field, i.e. to sin(g • x). 
Thus the Fourier transform of 5<p a picks up the components at ±<f. Finite size 
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of the domain spreads out the distribution, but if the domain has the size of, 
say, 3 wavelengths 27r/|g|, the half-width of peaks is only ~ 20% of \q\. Of 
course, larger domains lead to sharper peaks. If multiple domains contribute, 
then the signal should be folded with the momentum distribution of the do- 
main and averaged over directions. In Eq. ([]]) we have used the neutral classical 
pion field. This case leads to production of neutral pion quanta. However, we 
can rotate ansatz ([]]) in isospin, which leads to a degenerate solution. Such a 
configuration will decay into jets of charged pions. In addition, the quanta of 
the a field decay subsequently into pairs on neutral as well as charged pions. 
A more detailed calculation of the spectrum of pions requires additional as- 
sumptions for the dynamics of the phase transitions and is left for a further 
study. Also, the analysis of the formation and stability of domains is outside 
of the scope of this paper. 
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